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PREFACE 


This  report  develops  a  set  of  nonlinear  partial  differential  equations  which  describe  the  coupled  transfer 
of  energy  and  mass  through  hygroscopic  porous  media,  particularly  with  reference  to  textile  materials 
in  woven,  nonwoven,  and  laminated  forms.  The  report  focuses  on  the  derivation  of  the  equadons  and 
their  presentation  in  a  form  suitable  for  numerical  solution. 

This  work  was  undertaken  during  the  preliminary  preparation  of  a  doctoral  dissertation  proposal  to  the 
College  of  Engineering  at  the  University  of  Massachusetts  Lowell.  Dr.  Majid  Charmchi,  the  author's 
advisor  in  the  Department  of  Mechanical  Engineering,  provided  guidance,  suggestions,  and  encouragement 
during  the  derivation  of  the  equations  given  in  this  report.  The  two  technical  reviewers  at  the  U.S .  Army 
Natick  Research,  Development  and  Engineering  Center  were  Barry  Decristofano  of  the  Science  and 
Advanced  Technology  Directorate,  and  Gary  Proulx  of  the  Survivability  Directorate.  These  two 
reviewers  provided  especially  valuable  suggestions  and  comments,  and  pointed  out  several  errors  in  the 
draft  version  of  the  report. 


SYMBOLS  AND  ABBREVIATIONS 


Roman  Letters 
A  area  [m2] 

aop  surface  the  interface  per  unit  volume  [nr1] 

Ajt)  material  surface  [m2] 
c  constant  pressure  heat  capacity  [J/kg-°K] 

d  mass  fraction  weighted  average  constant  pressure  heat  capacity  [J/kg-°K] 

D  gas  phase  molecular  diffusivity  [nWsec] 

Da  diffusion  coefficient  of  water  vapor  in  air  [m2/sec] 

Dtff  effective  gas  phase  diffusivity  [m2/sec] 
g  gravity  vector  [m/sec2] 

h  enthalpy  per  unit  mass  [J/kg] 

h°  reference  enthalpy  [J/kg] 

h;  partial  mass  enthalpy  for  the  ith  species  [J/kg] 

heat  transfer  coefficient  for  the  a~P  interface  [J/sec-n^^K] 
enthalpy  of  vaporization  per  unit  mass  [J/kg] 
k  thermal  conductivity  [J/sec-m-°KJ 

K  [N/m1! 

d(Pc)/%T)  [N/m’-KJ 
K  permeability  coefficient  [m2] 

Kp  liquid  phase  permeability  tensor  [m2/sec] 

( msi )  mass  rate  of  desorption  from  solid  phase  to  liquid  phase  per  unit  volume  [kg/sec-m3] 

K/)  =  ^ l a*  Po  (p°  "  ^2 )  • 

(mw)  mass  rate  of  desorption  from  solid  phase  to  vapor  phase  per  unit  volume  [kg/sec-m3] 

(mh.)  mass  rate  of  evaporation  per  unit  volume  [kg/sec-m3] 

n  outwardly  directed  unit  normal 

p  pressure  [N/m2 } 

Pc  Py'Pf,'  capillary  pressure  [N/m2] 

p0  reference  pressure  [N/m2] 

p,°  reference  vapor  pressure  for  component  1  [N/m2] 

Q  volumetric  flow  rate  [m3/sec] 

Qtl  enthalpy  of  desorption  from  solid  phase  per  unit  mass  [J/kg] 
q  heat  flux  vector  [J/sec-m2] 


ix 


SYMBOLS  ANDABBREVIATIONS  (continued) 


f  position  vector  [m] 

r  characteristic  length  of  a  porous  media  [m] 

R.  gas  constant  for  the  ith  species  [N-m/kg-°K] 

T  temperature  [°K] 

T0  reference  temperature  [°K] 

T°  reference  temperature  [°K] 

T  total  stress  tensor  [N/m2] 

t  time  [sec] 

Uj  diffusion  velocity  of  the  ith  species  [m/sec] 
v  mass  average  velocity  [m/sec] 

V;  velocity  of  the  ith  species  [m/sec] 

VJt)  volume  of  the  solid  phase  contained  within  the  averaging  volume  [m3] 
VJ  t)  volume  of  the  liquid  phase  contained  within  the  averaging  volume  [m3] 
V  (t)  volume  of  the  gas  phase  contained  within  the  averaging  volume  [m3] 

•v  averaging  volume  [m3] 

VJt)  material  volume  [m3] 

h>  velocity  of  the  P~y  interface  [m/sec] 

u>j  velocity  of  the  a-y  interface  [m/sec] 

vv2  velocity  of  the  a~P  interface  [m/sec] 

Greek  Letters 

tg(t)  Vo/V,  volume  fraction  of  the  solid  phase 
E p(r)  Vp/'K  volume  fraction  of  the  liquid  phase 
ty(t)  \JV,  volume  fraction  of  the  gas  phase 

^  thermal  dispersion  vector  [J/sec-m3] 

t,  a  function  of  the  topology  of  the  liquid  phase 

|i  shear  coefficient  of  viscosity  [N-sec/m2] 

p  density  [kg/m3] 

p.  density  of  the  ith  species  [kg/m3] 

x  viscous  stress  tensor  [N/m3] 

x  toruosity  factor 

<l>  rate  of  heat  generation  [J/sec-m3> 

$  relive  humidity 

\  unit  tangent  vector 
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SYMBOLS  ANDABBREVIATIONS  (continued) 


Subscripts 

i  designates  the  ith  species  in  the  gas  phase 
/  liquid 

L  liquid 

s  solid 

5  solid 

a  designates  a  property  of  the  solid  phase 

P  designates  a  property  of  the  liquid  phase 

Y  designates  a  property  of  the  gas  phase 

ap  designates  a  property  of  the  a-p  interface 

cry  designates  a  property  of  the  o-y  interface 

Py  designates  a  property  of  the  P~y  interface 

Mathematical  Symbols 


didt 

DIDt 

d/dt 

<V> 


total  time  derivative 
material  time  derivative 
partial  time  derivative 

spatial  average  of  a  function  \| r  which  is  defined  everywhere  in  space 
phase  average  of  a  function  which  represents  a  property  of  the  P  phase 


intrinsic  phase  average  of  a  function  \j/p  which  represents  a  property  of  the  p  phase 
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Governing  Equations  for  Multiphase  Heat  and  Mass  Transfer  in  Hygroscopic  Porous  Media 

with  Applications  to  Clothing  Materials 


1.  Introduction 

The  purpose  of  this  report  is  to  develop  a  comprehensive  set  of  governing  equations  which  describe 
the  coupled  transfer  of  energy  and  mass  through  hygroscopic  porous  media.  The  report  focuses  on 
the  derivation  of  the  equations  and  ther  presentation  in  a  form  suitable  for  numerical  solution.  The 
various  steps  of  the  derivation  are  presented  in  sufficient  detail  so  that  the  origin  of  each  term  is 
clear.  The  level  of  detail  presented  should  make  it  easier  to  modify  each  equation  according  to  the 
requirements  of  a  particular  problem  or  material. 

The  basis  for  the  set  of  governing  equations  is  Whitaker's  comprehensive  theory  for  mass  and 
energy  transport  through  porous  media1.  These  equations  are  also  applicable  to  mass  and  energy 
transport  through  textile  materials  if  some  modifications  are  made.  Whitaker  modeled  the  solid 
portion  of  the  solid  matrix  as  a  rigid  inert  material  which  only  participates  in  the  transport  process 
through  its  thermal  properties.  In  hygroscopic  textile  materials  the  diffusion  of  water  into  the  solid 
is  a  significant  part  of  the  total  transport  process.  The  inclusion  of  the  extra  transport  terms  into  and 
out  of  the  solid  matrix  necessitate  extensive  modifications  of  Whitaker's  original  derivations. 

The  structure  of  this  report  follows  Whitaker's  derivations  as  closely  as  possible.  Many  references 
are  made  to  his  original  derivation.  Where  possible,  the  nomenclature  and  symbols  are  identical  to 
Whitaker's  original  derivation  to  facilitate  cross-referencing  between  this  modified  set  of  equations 
and  Whitaker's  original  set  of  equations. 
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2.  Mass  and  Energy  Transport  Equations 


The  hygroscopic  porous  media  is  modeled  as  shown  below. 


Figure  1.  Three  phases  present  in  hygroscopic  porous  media. 

A  typical  porous  hygroscopic  textile  material  may  be  described  as  a  mixture  of  a  solid  phase,  a  liquid 
phase,  and  a  gaseous  phase.  The  solid  a  phase  consists  of  the  solid  material  (usually  a  polymer  e.g.  wool 
or  cotton)  plus  any  bound  water  absorbed  in  the  solid  polymer  matrix.  The  solid  phase  is  thus  a  mixture 
of  the  solid  and  the  liquid.  This  definition  of  the  solid  phase  means  that  the  density  is  dependent  on  the 
amount  of  water  contained  in  the  solid  phase. 

The  liquid  (3  phase  consists  of  the  free  liquid  water  which  may  be  present  within  the  structure  of  the 
porous  solid.  This  would  also  include  water  which  is  contained  within  the  pore  spaces  of  the  solid  but 
is  not  sorbed  into  the  polymeric  matrix.  This  liquid  P  phase  is  a  pure  component,  and  we  will  be  able 
to  assume  a  constant  density  for  it. 

The  gaseous  y  phase  consists  of  the  vapor  component  of  the  liquid  (water  vapor)  plus  the  inert  air 
component.  Since  it  is  a  mixture  of  water  vapor  and  air,  its  density  will  not  be  constant,  but  will  be  a 
function  of  temperature,  concentration,  etc. 


The  appropriate  general  transport  and  conservation  equations  to  be  used  are: 


continuity  equation: 

|£+v-(pii)=o 


linear  momentum: 
Dv  -  r,  T 

p—  =pS+V‘T 


thermal  energy  equation: 

Dh  _  Dp 

p —  =  -V-5+— +Vv:  T+O 
v  Dt  H  Dt 

In  keeping  with  Whitaker's  derivation,  we  will  neglect  the  viscous  stress  tensor 


(2.1) 


(2.2) 


(2.3) 


Point  Equations 
a-Phase  -  Solid 


The  solid  a  phase  is  made  up  of  the  true  dry  solid  (polymer)  plus  any  of  the  liquid  phase  or  the 
vapor  component  of  the  gas  phase  which  has  dissolved  into  it  or  adsorbed  onto  its  surface.  This 
may  also  result  in  a  volume  change  for  the  solid  Q  phase  (swelling).  The  solid  is  now  a  mixture  of 
the  true  dry  solid  plus  the  liquid,  so  we  now  must  account  for  the  two  components. 


Since  swelling  is  possible,  which  results  in  a  small  velocity  of  the  solid  portion  due  to  its 
displacement,  we  must  account  for  the  velocity  of  the  solid  by  using  the  continuity  equation: 

3jk+V-(l :p.f,)-0  (2.4) 

and  for  the  two  components  of  liquid  (1)  +  solid  component  (2) ,  the  species  continuity  equation  is: 


^■+V,(p;v;)  =  °  ,  i  =  1,2,  ... 


dt 


(2.5) 


The  a  phase  density  is  not  constant,  since  it  includes  the  density  of  the  solid  plus  the  density  of  the 
liquid  contained  within  the  solid.  The  species  densities  are  likewise  not  constant,  since  the  species 
density  is  calculated  on  the  basis  of  the  total  phase  volume.  For  the  two  species: 


P  = 


=  Pi  +  p2 


(2.6) 


We  will  assume  that  the  dry  density  of  the  solid,  and  the  density  of  the  liquid  are  constant,  and  will 
call  them  p5  and  pL . 
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For  the  solid  phase,  we  can  divide  the  solid  phase  volume  into  the  fraction  taken  up  by  the  liquid, 
and  the  fraction  taken  up  by  the  solid: 


Volume  of  Liquid 
e°L  Total  a  Phase  Volume 


(2.7) 


The  relation  between  the  species  densities  and  the  solid  and  liquid  densities  is: 
P  =  eoLPz. +  (1  ~  tei)Ps  =  Pi  +  P2 
Pi  =  eglPl 
P2  =  (1~^JPs 


(2.8) 

(2.9) 

(2.10) 


The  density  and  velocity  of  the  mixture,  in  terms  of  the  species  densities,  is  given  as: 

Po=Pi+P2 

PoV0  =  p1Vi  +  p2v2 


(2.11) 

(2.12) 


or 

Po  =  &olPl  +  (1-£cl)Ps 
P  0*0  =  EolPl^I  +  (1  “  Eoi,)piSV2 


(2.13) 

(2.14) 


The  species  velocity  is  written  in  terms  of  the  mass  average  velocity  and  the  diffusion  velocity  as: 
vi=Va+ui  (2.15) 

and  the  continuity  equation  becomes: 


^+v.(p,.p0)=-v-(p,a,)  ,1=1,2, 3,... 


(2.16) 


The  diffusion  flux  may  be  written  in  terms  of  a  diffusion  coefficient  as: 

PA  =  ~P  0®ovf;r 

IPoJ 

and  the  continuity  equation  may  be  written  as: 

^+V.(p;v0)  =  V-|p0D0V^|  ,  i  =  l,  2,  3,  ... 


(2.17) 


(2.18) 


For  the  purposes  of  comparison  to  other  models  of  heat  and  mass  transfer  through  porous  materials, 
it  will  be  convenient  later  on  to  write  these  equations  in  terms  of  concentrations  of  water 
(component  1)  in  the  solid  (component  2). 

We  define  the  concentration  of  water  in  the  solid  ( Cs)  as: 
q  Mass  of  water  _  n\  _  p] 

s  Mass  of  the  solid  phase  ml  +  m2  pc  (2-19) 
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If  we  only  want  to  consider  the  continuity  equation  for  the  liquid  phase,  since  it's  really  the  only 
material  moving  into  or  out  of  the  solid  phase,  we  can  just  use  the  continuity  equation  for  the  liquid, 
which  is: 


3pi 

dt 


+V-(pi  v0)  =  V- 


\Po 


(2.20) 


Now  depending  on  how  we  want  to  treat  the  solid  velocity,  we  can  rewrite  this  a  couple  of  ways.  If 
we  say  that  the  solid  velocity  is  included,  then  in  terms  of  the  true  liquid  density,  where  the  species 
density  is  given  by: 


Pi  -  eojlPl  and  Pc  =  eoiPi  +  (1~  £<jl)Ps  > 


(2.21) 


the  continuity  equation  can  be  rewritten  as: 

"3er 
P  L 


or 


3ej 

dt 


3 1 


-+V-(eaLva)  = 


?+V-(£5lv(y)  =V.{pc%V(Cs)} 


v  PlJ  P  l 


(2.22) 


(2.23) 


If  we  can  neglect  the  solid  velocity,  the  continuity  equation  becomes: 
d£cL  _  f 


dt 


1  - ^  JV .[e0^L0V(Q)]  +  &{V .[®toV(Cj]} 
\  PlJ  P  l 


(2.24) 


We  must  also  include  the  momentum  balance: 

Pa^-  =  Paf  +  V-Ta  =>  Poj^  +  ^o -V)v0|  =  p0|  +  V-T0  (2.25) 


According  to  Jomaa  and  Puigali2,  we  may  also  write  the  linear  momentum  equation  as: 

Po'^  =  Pol+V'T°  (2,26) 
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There  are  a  couple  of  ways  to  address  the  mass  average  solid  phase  velocity.  If  we  assume  that  the 
total  thickness  of  the  materials  we  are  trying  to  model  does  not  change,  then  total  volume  remains 
constant,  and  the  change  in  volume  of  the  solid  is  directly  related  either  to  the  change  in  volume  of 
the  liquid  phase,  or  the  change  in  volume  of  the  gas  phase.  Another  approach  is  to  let  the  total 
volume  of  the  material  change  with  time.  As  the  material  dries  out,  and  the  tc*al  mass  changes,  the 
thickness  of  the  material  will  decrease  with  time  proportional  to  the  water  loss  which  takes  place. 
Allowing  the  thickness  of  the  material  to  change  with  time  would  result  in  a  solid  phase  velocity, 
which  we  could  relate  to  the  total  material  shrinkage.  The  two  situations  are  illustrated  in  Figure  2 
for  a  matrix  of  solid  fibers  undergoing  shrinkage  due  to  water  loss. 


=0 


Case  1 

Solid  liber  shrinkage 
results  in  bulk  thickness 
reduction  and  nonzero 
mass  average  solid  velocity. 


Case  2 

Total  bulk  thickness  and  volume 
do  not  change;  shrinkage  of 
solid  fiber  portion  due  to 
water  loss  does  not  result  in  a 
mass  average  solid  velocity. 


Figure  2.  Two  methods  of  accounting  for  shrinkage/swelling  due  to  water  uptake  by  a  porous  solid. 

We  will  assume  that  the  shrinkage  behavior  is  like  the  first  case  shown.  This  means  that  we  must 
include  the  mass  average  velocity  in  the  derivations,  and  that  the  total  material  volume  (or  thickness 
in  one  dimension)  will  no  longer  remain  constant. 

Jomaa  and  Puiggali  also  give  an  equation  for  the  solid  velocity,  in  terms  of  the  intrinsic  phase 
average  (discussed  later)  as: 

Ma  =  {p)'%-rJo|'<pa>^  <2-27> 

where  £  is  the  generalized  space  coordinate,  with  the  origin  at  the  center  of  symmetry,  and  n 
depends  on  the  geometry  (n=l -plane, n=2~cylinder,  n=3~  sphere)  according  to  the  paper  by 
Crapiste  et.  al.3 
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The  thermal  energy  equation  is: 


(2.28) 


Some  simplifying  assumptions  can  be  made  at  this  point  by  neglecting  several  effects.  We'll  start 
by  dropping  the  reversible  and  irreversible  work  terms  in  the  thermal  energy  equation,  along  with 
the  source  term,  and  expand  the  material  derivative: 


(2.29) 


It  will  be  assumed  that  enthalpy  is  independent  of  pressure,  and  is  only  a  function  of  temperature, 
and  that  heat  capacity  is  constant  for  all  the  phases. 

We  can  replace  the  enthalpy  by: 
h  =  cpT  +  constant  ,  in  the  c,  p,  and  y  phases. 


We  can  now  rewrite  the  thermal  energy  equation  as: 

3  M„r< 


Pe¬ 


ar 


•+p, 


r  vo  -V  (c>)c 


Ta  +  constant 


(2.30) 


(2.31) 


(2.32) 


Pa(^)0‘^  +  Po(cp)0  V^a  Mcp)0{*^  + 


We  may  apply  Fourier's  law  to  obtain: 
po(^)0{^+foVro}  =  *„V:T0 

or,  for  a  multi -component  mixture: 

Pa(e,) Il^+VVJiW’j'-V.  Xp jBjhj 
V  0  J  U=1  / 

j=Np  . 

where  (cp)a  = 

j= 1  Pa  1 

and  the  partial  mass  heat  capacity  and  enthalpies  (cp).  ,  hj  are  given  by  the  partial  molar 
enthalpy  and  the  partial  molar  heat  capacity  divided  by  the  molecular  weight  of  that  component. 


f  i=N 


(2.33) 
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P  Phase  ••  Liquid 


The  continuity  equation  for  the  liquid  phase  is: 


if+v  to>)=° 


(2.34) 


For  the  thermal  energy  equation,  as  we  did  before,  we  neglect  compressional  work  and  viscous 
dissipation: 

-2£  =  Vvp:ip  =  <I>p=0  (2. 

which  reduces  the  thermal  energy  equation  to: 


(2.35) 


(dh*  \ 

Pp  =“v^p 

\  dt 


(2.36) 


If  we  assume  enthalpy  only  depends  on  temperature  and  specific  heat,  as  we  did  for  the  solid,  we 
may  write  the  thermal  energy  equation  for  the  liquid  phase  as: 


pph)p(^+i,PVJi.)=^ 


(2.37) 


The  liquid  momentum  equation  will  be  discussed  later  in  terms  of  a  permeability  coefficient  which 
depends  on  the  level  of  liquid  saturation  in  the  porous  solid. 
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y  Phase  -  Gas 

The  gas  phase  is  made  up  of  the  vapor  form  of  the  liquid  (5  phase,  and  an  inert  component  (air).  We 
do  not  need  to  modify  any  of  the  assumptions  made  by  Whitaker  for  this  phase,  so  we  may  simply 
write  down  the  equations  given  by  Whitaker1: 

continuity  equation: 

^+V(pTvy)  =  0  (2.38) 

and  for  the  two  components  of  vapor  (1)  +  inert  component  (2) ,  the  species  continuity  equation  is: 

%i-+V.(o,.v(.)  =  0  ,  i  =  1,2,  ...  (2.39) 

dt 


The  density  and  velocity  of  the  mixture  is  given  as: 
Py  =  Pi +  P2 
PY?Y=P1?1  +  P2?2 


(2.40) 

(2.41) 


The  species  velocity  is  written  is  terms  of  the  mass  average  velocity  and  the  diffusion  velocity  as: 
Vi^Vy+tt,  (2.42) 

and  the  continuity  equation  becomes: 


V-(p;vY)  =  -V-(p;w() 


,  i  =  1,  2,  3, 


(2.43) 


The  diffusion  flux  may  be  written  in  terms  of  a  diffusion  coefficient  as: 


P  fit  =  ~Py®V 


Pi_ 

Py 


(2.44) 


and  the  continuity  equation  may  be  written  as: 


3a 

dt 


L+V.(p,vy)  =  V. 


' 

(n  ^1 

Py2*7 

£l 

IPyJ. 

,  i  =  1,  2,  3,  ... 


(2.45) 


It  is  possible  that  we  can  neglect  the  change  in  gas  density  with  time,  or  at  least  the  change  in  the 
density  of  the  inert  component,  and  only  consider  the  continuity  equation  for  the  vapor  component 
of  the  gas  phase  (component  1): 

f  W 

Pl 


|L+V.(PlfT)  =  V- 


PyPV 


(2.46) 


and  if  we  have  no  gas  phase  convection,  with  the  gac  phase  stagnant  in  the  pore  spaces,  the 
continuity  equation  becomes: 


3a=vJ 

dt 


py2?V 


\ 

Pl 

VPT  Jj 


(2.47) 
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The  thermal  energy  equation  is  given  as: 


•VTv 


-ky  v2r-v 


fi=N  _  ^ 

£pM- 

V  1=1 


where 


(2.48) 


and  the  partial  mass  heat  capacities  and  enthalpies  ( cp ).  ,  ht  are  again  given  by  the  partial  molar 
enthalpy  and  the  partial  molar  heat  capacity  divided  by  the  molecular  weight  of  that  component. 


10 


Boundary  Conditions 

Whitaker  next  derives  the  boundary  conditions  for  each  phase  interface.  This  section  of  the 
original  derivation  must  be  extensively  modified  since  we  no  longer  have  a  rigid  solid  phase  with 
zero  velocity.  We  will  no  longer  have  a  simple  set  of  boundary  conditions  for  the  solid-liquid  and 
solid-vapor  interface.  The  conventions  and  nomenclature  for  the  phase  interface  boundary 
conditions  are  given  below  in  Figure  3,  which  follows  Whitaker’s  approach  as  closely  as  possible. 


Figure  3.  Material  volume  containing  a  phase  interface,  with  velocities  and  unit  normals  indicated. 
Only  two  phases  (solid  and  gas)  shown. 


Liquid-Gas  Boundary  Conditions 

Whitaker  gives  the  appropriate  boundary  conditions  for  the  liquid-gas  interface  as: 


»=N 


(vp  -  )  •  j «pY  +  Pyhy  (vY  -  M>)  •  =  -  %  ■ %  +  +  X  PM  •  «yfi  ’ 

(2.49) 

pp(iTp  -w)-^  +pY(?Y  -*)■«£  =0 

continuous  tangent  components  to  the  phase  interface  h 

(2.50) 

•C?' 

II 

-r1 

(2.51) 

species  jump  condition  given  by: 

P/(v/-w>)-nYp+pp(vp-Vp)-^Y=0  ,  i  =  l 

Pl(v,  -vv)-HYp  =  0  .  i  =  2,3,  ... 

(2.52) 

(2.53) 
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Solid-Liquid  Boundary  Conditions 

The  boundary  conditions  for  the  solid-liquid  interface  are  identical  except  that  the  phase  interface 
velocity  is  given  by  w,  . 


PcAt  0 Vo  -  *2  )  ■ * ' "op  +  Pp^p  ( Vp  -  A )  •  «pc  " 

M 

"op 

P  o  ( -  ^’2 )  •  +  Pp  ( ~  ^2 ) '  =  0 

continuous  tangent  components  to  the  phase  interface  X.: 

Vo  ‘  ^o(3  =  Vp  ^po 

species  jump  condition  given  by: 

Pj  ( ”  *2 ) ' 1  «po  +  Po  ( ~  *2 ) ' ‘  "op  *  0  ’  j  =  1 

P;(*7-*2)'V  =  °  ’  J  =  2’3’  - 


(2.54) 

(2.55) 

(2.56) 


(2.58) 


Solid-Gas  Boundary  Conditions 


The  boundary  conditions  for  the  solid-liquid  interface  are  modified  even  more  because  we  have  a 
phase  interface  between  two  multi-component  phases.  The  phase  interface  velocity  is  given  by  wt . 


Po^o (^o  “  ^1 ) ’ ^oy  ■*" Py\  ( Vy  ~  *^1 ) Hya ~  j 


$a  +  ^,PyffA 

;=1 


”oy 


i-N 

i=l 


(2.59) 


Pa^o-W'OAry  +  pY^y-^l)  'V  =0 

continuous  tangent  components  to  the  phase  interface  X: 

Vo  ‘  ^oy  =  Vy  ■  ^yo 

species  jump  condition  given  by: 

P;(^-vv1)-«oy+P.(v,  -M>i)  Ao  =°  -i  =  1  ’J  =  1 
p;(v;  -vPi)-/iOY=0  ,  j  =  2,3,  ... 

Pi(^-^i)*V  =0  .  i  =  2-3-  - 


(2.60) 


(2.61) 


(2.62) 

(2.63) 

(2.64) 
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Volume  Averaged  Equations 


Whitaker  uses  the  volume-averaging  approach  outlined  by  Slattery4  so  that  many  of  the  complicated 
phenomena  going  on  due  to  the  geometry  of  the  porous  material  are  simplified.  He  defines  three 
averages: 


Spatial  Average:  Average  of  some  function  everywhere  in  the  volume. 

Phase  Average:  Average  of  some  quantity  associated  solely  with  each  phase. 

Intrinsic  Phase  Average: 

ft 

O 


We  can  also  define  volume  fractions  for  the  three  phases  as: 


eo(0  = 


Ml  .  ,rt.M 


Ep(0  =  ' 


ey(/)  =  - 


(2.65) 


(2.66) 


(2.67) 


(2.68) 


We  now  have  the  volume  and  volume  fraction  of  the  solid  changing  with  time,  which  is  not  the  case 
with  Whitaker's  original  derivation. 


We  are  going  to  say  that  the  total  volume  is  conserved,  or  that: 

^-Vo(0  +  Vp(0  +  VT(0 

The  volume  fractions  for  the  three  phases  are  related  by: 

eo(0+ep(f)  +  eY(0-l 

and  the  phase  average  and  the  intrinsic  phase  averages  are  related  as: 

s.W'=W 


(2.69) 


(2.70) 


(2.71) 
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Volume  Average  for  Liquid  (3  Phase 


We  will  first  look  at  the  volume  average  for  the  fi  phase.  It  will  be  complicated  because  of  the  three 
different  phase  interface  velocities  which  we  must  now  include  in  the  analysis,  whereas  Whitaker 
only  had  to  account  for  the  single  liquid-gas  interface  velocity. 


The  continuity  equation  for  the  liquid  phase  is: 


(2.72) 


We  will  integrate  over  the  time-dependent  liquid  volume  within  the  averaging  volume,  and  divide 
by  the  averaging  volume  to  obtain: 


We  may  take  the  first  term: 


-j 

<k'JvP(*)V 


and  apply  the  general  transport  theorem3: 
d  r  r  cty 

•'(D  or  -^(i) 


(2.73) 


(2.74) 

(2.75) 

(2.76) 


We  note  that  *F  = 

Bt 

and  using  the  modified  general  transport  theorem  we  obtain: 

(3p  ^ 

it  f =  7t  ?  JifrM “ ?  L* PP* ‘ V ^  t  h  PP*2  '  V44  (2.77) 

For  the  other  term: 

we  may  use  the  volume  averaging  theorem: 

NO  =  VN}+ ^  Vl ^M+^r 

to  rewrite  the  term  as: 

1  Km  V  ’  =  ( V '  (PPPP ))  =  V  ’ (pf^P ) +  ^  J^()  Pp*P '  + ^  1V(()  pPpP 1 V*4  (2.80) 


(2.78) 


(2.79) 
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Whitaker  says  that  we  may  also  rewrite  the  time  derivative  term  as: 


d_ 

dt 


d_ 

dt 


d_ 

dt 


This  allows  us  to  rewrite  the  continuity  equation  for  the  liquid  phase  as: 
|(pp)  +  v.{ppvp)  +  ij^p^vp -4^  +  i]^Pp(vp =  0 

We  may  assume  that  the  density  in  the  liquid  is  constant,  so  that: 

(ppVp)  =  Pp(vp) 

(pp)  =  £pPp 


(2.81) 


(2.82) 


(2.83) 

(2.84) 


The  liquid  velocity  vector  may  be  used  to  calculate  volumetric  flow  rates.  Whitaker  gives  the 
example  of  the  flow  rate  of  the  liquid  phase  past  a  surface  area  as: 

Qf=jA(vf)-ndA 


(2.85) 


The  two  constant-density  liquid  relations  given  above  allow  the  liquid  phase  continuity  equation  to 
be  rewritten  as: 


(2.86) 


The  thermal  energy  equation  for  the  liquid  phase  was  given  previously  as: 


Pp 


=  -V<7, 


% 


Whitaker  notes  that  this  may  be  rewritten  by  adding  the  term 

side  to  get: 

d 

jt  (ppfy ) + v ■  (pp  Vp  ) = v  -  % 


%+v-(prt) 


(2.87) 


to  the  left-hand 


(2.88) 
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We  may  follow  the  same  procedure  used  previously  to  obtain  a  volume  averaged  form  of  the 
thermal  energy  equation  where  we  use  the  general  transport  theorem  for  the  first  term  and  the 
averaging  theorem  for  the  second  and  third  terms  to  get: 

^(pPV)+v-(pe/rt)+  vL/p'fcte 

=  (2.89) 

All  we  did  here  was  add  an  additional  term  to  Whitaker's  equations  due  to  the  solid-liquid  interface 
velocity. 


Whitaker  uses  the  relation  for  the  enthalpy  of  the  liquid  phase: 

fy=V+Mp(rP“rp‘)  (2-9°) 

and  goes  through  several  steps,  accounting  for  the  deviation  and  dispersion  terms  from  the  average 
properties  (marked  with  a  tilde),  to  write  an  expression  for  the  two  terms: 


+ V  ■(pp^%)  =  epp^c,  )p 


dWL+ p 

B,  +P» 


v+^)p((rp)p-3p-)|^+v.(pp) 


+Pp(cp)p{^p}‘^(2Ji)  +Pp(cp)p^‘(7p^p) 


(2.91) 


We  can  recognize  that  the  term: 


3e 


lf+vW 


is  contained  in  the  liquid  phase  continuity  equation: 

^  +  V-(9p)]+^J^j^-vv)-nprd4+lj^(^-w2).^(i4  =  0 


at 

\ 

so  that: 


3 

The  expression  for  the  two  terms  '^(pp^p)  +  ^,{pp^1p^p)  may  be  written  as: 

^■(pp^p) + v  •  (pp  Vp )  =  Eppp  (cp  )p  ■  +Pp  ( cp  )p  (  Vp  )  •  V(rp  f  +  pp  ( cp  )p  V  •  (fpfp  ) 


(2.92) 


(2.93) 


V  +  (cp)p((7p)P-7P#)  +  (2.94) 
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We  may  now  substitute  back  into  the  thermal  energy  equation  for  the  liquid  phase: 

bIt  \P 

eppp(cf )p  '[''^'-+Pfi(cpk{''ls)'V{Tti) 

L*  Pp  1 ^ )p  (rp  ”  X*fc  -  *) ' 1 Hto‘tA 

+?J 

-  v)(V "«»)■ 'V“ 

=  -V •{?„) +{®p>  < %  ■  %dA ~^'V«  (2.95) 

Whitaker  now  uses  Gray's  definitions  of  the  point  functions  for  the  phase  properties6  as: 

rp=(rp)p+fp  (2.96) 

This  allows  the  liquid  phase  thermal  energy  equation  to  be  written  as: 

EpPp  (c,.)p  +pe  {% )  ‘  v{rp  )P + Pp  <cp  v '  ) 

=  (2.97) 

Whitaker  now  rewrites  the  heat  flux  term  using  Fourier's  law  (fy  =  -ApVTp)  and  the 

averaging  theorem  to  write: 

<?b)  -  -*p<V3i>)=-*p  V(7P>+5;Lfta  (2.98) 
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We  may  also  substitute  the  intrinsic  phase  average  temperature 

for  the  phase  average  temperature  (t^  to  obtain  an  expression  for  the  heat  flux  vector: 


The  thermal  energy  equation  for  the  liquid  phase  may  now  be  written  as: 


c(lP[lfcp  k  “jT”  +P»  <Cf  ) '  V{rP  }P  +  P»  (cP  J0  7 '  {tyt) 


=v- w 


vfeafcfV-f  r0S..<M+— f  T*RdA 
V  P'  P'  J  1,JV  p  *°  <V‘%  J 


-£lVlvvA-iL%-vM 


(2.99) 


(2.100) 
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Volume  Average  for  Gas  y  Phase 


The  gas  phase  continuity  equation  is  identical,  up  to  a  point,  to  the  continuity  equations  for  the  solid 
and  liquid  phases: 

3  1  1 

^(pr)+V-{pTv1>+-J^pT(vY-»)-^i4+-JiVPY(vT-Hi1).^<14  =  0  (2.101) 

For  the  liquid  and  solid  phases  we  could  assume  constant  density,  and  simplify  the  equation  further; 
we  can't  do  this  for  the  gas  phase  since  the  density  depends  on  the  temperature  and  the  pressure. 

Whitaker  uses  the  Gray's  expressions  for  the  point  functions  again,  along  with  the  definition  of  the 
intrinsic  phase  average  to  rewrite  the  gas  phase  continuity  equation  as: 

|(eJ(pT)7)+V-({pr)1'(vT))+V-{pYaY> 

+  (2.102) 

Whitaker  then  assumes  we  can  neglect  terms  with  products  of  the  dispersion  or  deviations,  so  we 
can  drop  that  term  to  write  the  gas  phase  continuity  equation  as: 

|(EY{pT)T)+v{(pr)1'(vT))+^JjVPT(j’T-«>)-^+ij^p,(i>T-».1)  BIad4  =  0  (2.103) 

Since  the  gas  is  a  multi-component  mixture  we  must  also  go  through  the  species  continuity 
equation: 

3  1  1  . 

0^<p«> +  V ■  <P| V,-) + ~  J  P; (v,  -  w) •  rLfidA+—j^ p,(v,  -  vPj )  fL^dA  =  0  i  =  1,2,  ...  (2.104) 

Whitaker's  derivations  may  be  used  directly,  and  we  can  write  the  final  form  of  the  gas  phase 
species  continuity  equation  as: 

(e Y <P/ > Y )  +  V •  (<P/ > Y ( v( >)  +  “ Pi (v,  - w) ^ p,- ( v(- -  w, ) - ^4  (2.i05) 

If  we  neglect  the  deviation  terms,  and  if  we  also  consider  only  the  species  continuity  equation  for 
the  vapor  component  (component  1),  we  can  rewrite  the  gas  phase  continuity  equation  as: 

f  V 

|K(Pl)’')+V  «P,)T(vT))+^/^Pl(vi-«1)  Slf)*l  =  V  '(PT)T®V  ■  (2.106) 
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The  corresponding  thermal  energy  equation  for  the  gas  phase  may  also  be  written  as: 


f  i=N 

i=N 


pfl>  -w 


5> i(cp)ify{vi-w)  n^dA 


+'!p\Ayg  Sp«' Shitty  -  ) •  fi^dA 

i=i  1=1 


=V.^ 

If 


V(er{rr)  )  +  v 7r'iKJlM  +  ,,  It,  7yr‘ifi‘14 


'I'' 


(2.107) 


20 


Volume  Average  for  Solid  c  Phase 

The  volume  averaging  procedure  for  the  liquid  phase  was  made  general  enough  so  that  the  same 
equations  also  apply  for  the  solid  phase.  The  only  differences  are  that  now  the  phase  interf  ace 
velocities  are  w2  for  the  solid-liquid  interface,  and  wx  for  the  solid-gas  interface.  We  also  need  to 
account  for  the  species  continuity  equations.  Since  the  two  components  (liquid  and  the  solid)  are 
assumed  to  have  a  constant  density,  we  will  not  run  into  the  same  complications  we  did  with  the  gas 
phase  continuity  equation.  The  appropriate  subscripts  for  the  solid  phase  also  need  to  be  added  to 
the  equations. 

We  cannot  assume  that  the  solid  phase  density  is  constant,  since  it  is  a  mixture  of  the  solid  and  the 
liquid  component.  However,  it  will  be  less  complicated  than  the  gas  phase  density  since  we  can 
assume  that  each  component's  density  is  constant. 


The  solid  phase  continuity  equation  is: 

f  .1-.  \ . 


xy  tvjuauuu  to. 

and  the  species  continuity  equation  is: 

•  i  -  -  -  1 


(2.108) 


ana  me  species  cuuuiiuuy  cijuauvm  i*. 

|((p> )  +  ^ ' (Pv^/ )  +  ^ ~  ^ +  “7  {$j  -w2)-napdA  =  0  j  =  l,  2,  ...  (2.109) 

We  can  follow  the  same  derivation  used  for  the  gas  phase  to  write  the  gas  phase  continuity  equation 
as: 

and  the  final  form  of  the  solid  phase  species  continuity  equation  is: 

r  r  /  yi  1 


(2.110) 


=  V- 


i' (p°>x  wJ 


j  =  1*2,  ... 


(2.111) 


If  we  want  to  just  follow  the  single  liquid  component  (component  1)  and  write  the  continuity 
equation  for  that  species,  we  may  write: 

|K<Pir)+v((p,r<v1))+^p,(n-»2)AS‘M+^j^Pi(iii-n)«a,‘M 

f  M 


=  V- 


(Po)°^V 


Pi 


,(Pc)°JJ 


(2.112) 
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Later  on,  we  may  also  want  to  assume  that  the  solid  velocity  is  zero,  so  we  could  rewrite  the  solid 
phase  continuity  equation  as: 


=v- 


pi 


v(Po)  J 


(2.113) 


The  corresponding  thermal  energy  equation  for  the  solid  phase  may  also  be  written  as: 

JE(P;><cp)y 

;=1 


IMJ 


3r 


j=N 


M 


1  .  W 


+  vSa^  2Pj(cp')M*J-*iYfi<rtdA 
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(2.114) 


This  completes  the  continuity  and  thermal  energy  volume  averaged  equations  for  all  three  phases. 
The  various  continuity  equations  are  given  in  several  forms,  depending  on  whether  we  want  to 
include  the  solid  velocity,  and  whether  we  just  want  to  use  the  continuity  equation  for  the  liquid 
component  only,  since  it  is  the  only  species  which  is  transferring  between  the  three  phases. 
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3.  Total  Thermal  Energy  Equation 


The  three  phases  are  assumed  to  be  in  local  thermal  equilibrium  so  that: 
(Ta)° ={Tyy 

(T^t^y +H(Tsif +ty(Tyf =(r„)“  .(Ttf =(r1)T 


(3.1) 

(3.2) 


We  can  now  write  the  total  thermal  energy  equation  by  adding  together  the  thermal  energy 
equations  for  each  phase,  and  using  the  local  thermal  equilibrium  relations  given  above.  This 
equation  is  identical  to  Whitaker  s,  except  for  the  addition  of  extra  terms  due  to  the  solid-gas  and 
solid-liquid  phase  interface  velocities,  which  are  no  longer  zero.  The  equation  is  also  written  so 
positive  flux  terms  imply  liquid  evaporating  into  the  gas  phase,  rather  than  condensing. 
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(3.4) 


Whitaker  defines  a  spatial  average  density: 


j=N 


i=N 


<P> = <=o  S(py)°  +eu(pp)  +eT  X(p,)r 

;=i  i=i 

and  a  mass  fraction  weighted  average  heat  capacity  by: 

e«;i(pyr(c?);+^PP(^)p+eY'f(pi)v(c,) 

c. — e! - - - <-■ 

'  <p) 

This  allows  the  first  term  in  the  thermal  energy  equation  to  be  written  as: 

[i-N 


(3.5) 


U=i  J 


' +  EpPp  [Cp  )p  +  eY  |  X  <P«  )(C/>  );  | 
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(3.6) 


We  must  now  consider  the  interphase  flux  terms  in  the  total  thermal  energy  equation.  In  Whitaker’: 
derivation,  he  only  had  one  inteiphase  flux  term  to  consider,  that  of  the  exchange  of  mass  between 
the  liquid  and  the  gas.  We  now  have  two  more  interphase  fluxes  to  consider:  that  between  the 
liquid  and  the  solid,  and  that  between  the  gas  and  the  solid. 

We  first  follow  his  derivation  for  the  liquid-gas  interface,  and  then  apply  it  to  the  other  two 
interfaces. 

The  jump  boundary  condition  for  the  liquid-gas  interface  was  shown  previously  to  be: 
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Ppfy  (vp  -  w)  ■  /ip7  +  Py/ly(vy  «  *)  ■  =  ~  ’  % 

and  this  may  be  rewritten  as: 
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The  jump  boundary  condition  for  the  solid-gas  interface  was  shown  previously  to  be: 

PoAo(^o-lS>l)-««Ty+pYAy(i,y-w1)ftyo  =  - 
and  this  may  be  rewritten  as: 
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The  jump  boundary  condition  for  the  solid-liquid  interface  was  shown  previously  to  be: 
P  A  (vc  - 1 vi>2 ) •  ■ «<jp  +  P  A  ( Vp  “  % )  •  «p<, =  ~ 

This  may  be  rewritten  as: 
j=N 

IPAfe  -*2)  Ap +pp%(vp  -  vv2)-«p<j = -(% 

;'= 1 


(3.11) 


Ap 
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Using  these  results,  we  may  write  the  interphase  flux  terms  in  the  total  thermal  energy  equation  as: 
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The  total  thermal  energy  equation  is  now  written  as 

Hi  im  s 

niZxlV/,  L  W,\, 

'P 


<p 

-If 

<V  J/Lo 


Z(cP);{p/;>  +  Pp(^)p{vp)+  X(cp).{p;V(> 

L  i=i 

i=N  r  _  1 

X  Py  [ hj  -  (CP  )j  ?c  J( ' ~  ^2 )  •  ^jp 

[+Pp[^-(cp)pfp](vp-vv2)-np0 

Pp[ty~(Cp)p%](fy“*),'fcy 


W) 


-if 


.if 

<77  J/ 


»=1 


X  p;  [ty  -  {CP  )j  ?c  ](vy  -  VV! )  •  «0Y 

i=N  r-  -  i 
h  X P*[^*»  ~{cp  );  Ty  J(V|  “  )  •  Hyo 


>dA 


*V- 


[v[(Mo+Vp+VyKr>] 

^{ka~^)l^\^T^dA 

+(^o  -  *y  )  ^7  Ty  n<r,dA 


(3.14) 


26 


Next,  we  can  begin  to  express  the  phase  interface  velocities  in  terms  of  enthalpies  of  vaporization, 
sorption,  and  desorption. 

The  enthalpies  for  each  phase  were  defined  previously  as: 


© 

f-45 

1 

+ 

© 

•sT 

ii 

l-sT* 

(3.15) 

*P=V+(cp)p(rP-¥) 

(3.16) 

*-V +(<>), ft-  Ty‘) 

(3.17) 

We  also  know  that  the  intrinsic  phase  average  temperatures,  temperature  dispersion,  and  overall 
average  temperatures  are  related  by: 


(t,)'  =(!•„)"  =<rT)T  ={r> 


We  may  use  these  relations  to  rewrite  the  integrands  inside  the  volume  integrals  on  the  left  hand 
side  of  the  total  thermal  energy  equation.  Whitaker  gives  the  result  for  the  liquid-gas  interface  as: 


(3.22) 


From  the  species  jump  conditions: 

P«  (Vf  -  vv)  •  +  pp(vp  -  w)  •  rtp-y  =  0  ,  i  =  1  (3-23) 

Pi(vi-w)  n^  =0  ,  i  =  2,3,  ...  (3.24) 

where  the  subscript  1  refers  to  the  component  (water)  which  is  actually  crossing  the  phase  boundary 
as  it  goes  from  a  liquid  to  a  vapor. 
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From  the  species  jump  conditions  we  may  also  write: 
P,  ( V!  -  W2 )  •  na p  =  -pp  (^P  -  H>2  )  •  "Per 


dA 


We  may  rewrite  the  integral  as: 

_lf  |[*P -(^)p(<r>-rP )]pp(% 
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-  -  (cp  )i  «7’>  -  ^  )]pi  -*H» 
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(3.25) 


(3.26) 


We  may  use  the  following  definitions: 


A Kv  (at temperature  (T))  =  {[*;, -^+(cp\((.T)-T-)-(cp)^T)-T^ 

K)  =  ^pp(fp--)^ 


to  rewrite  the  integral  as: 

,  Pfif-V  -  (C/>  )p  fp](F0  -  ^  -* 

V  **  ip.[*.-(tr),7r](f1--)  'V 

.  1=1 


\M=&hvap(mlv) 


(3.27) 

(3.28) 


(3.29) 


The  corresponding  terms  for  the  phase  interface  between  the  solid  and  the  liquid  are  identical, 
except  that  we  no  longer  use  the  quantity  A h  ,  but  instead  use  the  differential  enthalpy  of 
sorption7,  which  we  will  give  the  notation  Q , .  The  differential  heat  of  sorption  is  the  heat  evolved 
when  one  gram  of  water  is  absorbed  by  an  infinite  mass  of  the  solid,  when  that  solid  is  at  a 
particular  equilibrated  moisture  content.  This  is  very  similar  to  the  heat  of  solution  or  heat  of 
mixing  that  occurs  when  two  liquid  components  are  mixed.  For  textile  fibers  there  is  a  definite 
relationship  between  the  equilibrium  values  of  the  differential  heat  of  sorption  and  the  water  content 
of  the  fibers,  and  we  can  use  those  relationships  in  our  thermodynamic  equations  which  will  be 
discussed  in  a  later  section. 


28 


The  solid-liquid  interface  integral  term  is  thus  given  as: 
j=n  r_  i  ^4 
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From  the  species  jump  conditions  we  may  also  write: 

Pi  ft  “  *2 )  ■ ■  ■ «ap  =  “Pp  (%  - *i )  ‘  «po 

We  may  rewrite  the  integral  as: 
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Qx  (at  temperature  (T))  = 
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to  rewrite  the  original  integral  as: 
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U4  =  G,(m„)  (3-35) 
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For  the  gas-solid  interface,  the  heat  of  desorption  for  the  vapor  is  equal  to  the  energy  required  to 
desorb  the  liquid  plus  the  enthalpy  of  vaporization  required  to  evaporate  the  liquid: 

Q^Qt+AKap  (3.36) 


The  derivation  is  exactly  the  same  as  for  the  other  two  interfaces,  where  the  only  component 
crossing  the  phase  interface  is  component  1  (water)  and  we  may  write  the  integral  as: 


—  ((?/  &hvap  ){^(rv  ) 


(3.37) 


For  these  equations  { msi )  is  the  mass  flux  desorbing  from  the  solid  to  the  liquid  phase,  («!„)  is 

the  mass  flux  desorbing  from  the  solid  into  the  gas  phase,  and  (m^)  is  the  mass  flux  evaporating 
from  the  liquid  phase  to  the  gas  phase. 


The  total  thermal  energy  equation  now  becomes: 
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We  may  simplify  the  total  thermal  energy  equation  based  on  an  effective  thermal  conductivity,  and 
write  our  total  thermal  energy  equation  in  a  much  shorter  form  as: 


+Ahvap  (mlv)  +  Qi(msl)  +  (Q,  +  A h^ap  \m^) 

=  V-(^-V(r))  (3.39) 


The  effective  thermal  conductivity  can  be  expressed  in  a  variety  of  ways  as  described  by  Whitaker, 
depending  on  the  assumptions  you  choose  to  make  about  the  isotropy  of  the  porous  medium,  the 
importance  of  the  dispersion  terms,  etc.  The  effective  thermal  conductivity  is  also  an  appropriate 
place  to  include  radiative  heat  transfer,  and  one  could  add  an  apparent  radiative  component  of 
thermal  conductivity  to  the  effective  thermal  conductivity  to  account  for  radiation  heat  transfer. 
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4.  Thermodynamic  Relations 


The  gas  phase  is  assumed  to  be  ideal,  which  gives  the  intrinsic  phase  partial  pressures  of  the  gas 
phase  as: 


W  =(pifKi(n  i  =  u.  - 


(4.1) 


We  also  have  the  relations  for  the  gas  phase,  where  for  our  case  component  1  is  water,  and 
component  2  is  air: 


(Py)Y  ={Pi)Y+{P2>Y 
{Pyf  ={Pl)y  +(P2? 


We  must  also  connect  the  differential  heat  of  sorption,  Qr  with  the  concentration  of  water  in  the 
solid  phase.  An  example  of  a  general  form  for  Qr  ( in  Joules/gram)  can  be  expressed  as  a  function 
of  the  relative  humidity8: 

r  i  i  \  _  /_\Y 


Qt  (J/g)  =  195(1  -<(>) 
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-+ 
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(0.2 -'•<}>)  (1.05-4)) 


where  relative  humidity  4>  =  —  = 

Ps  Ps 


(4.4) 


<(>  Relative  Humidity 

Figure  4.  Generic  differential  heat  of  sorption  for  textile  fibers  (sorption  hysteresis  neglected). 

We  must  connect  the  differential  heat  of  sorption  with  the  actual  equilibrium  water  content  in  the 
solid  phase.  For  the  two  component  mixture  of  solid  (component  2)  plus  bound  water  (component 
1)  in  the  solid  phase,  the  density  of  the  solid  phase  is  given  by: 

(Po)°={Pi}0'+{P2)°  (4.5) 
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We  could  make  the  assumption  that  mass  transport  in  the  textile  fiber  portion  is  so  rapid  that  the 
fiber  is  always  in  eqnili’ '  ;,im  with  the  partial  pressure  of  the  gas  phase,  or  is  saturated  if  any  liquid 
phase  is  present.  This  wo.  rUminate  the  need  to  account  for  the  transport  through  the  solid  phase 
at  all.  There  are  a  variety  of  sorption  isotherm  relationships  we  could  use,  including  the 
experimentally-determined  relationships  for  a  specific  fiber  type,  but  a  convenient  one  is  given  by8: 


Regain  (R)  =  fy(0.55<)>) 


1  1 
(0.25  +  (J))  +  (1.25-<{)) 


(4-6) 


Rf  is  the  standard  textile  measurement  of  grams  of  water  absorbed  per  100  grams  of  fiber,  measured 
at  65%  relative  humidity.  We  may  rewrite  this  in  terms  of  the  intrinsic  phase  averages  for  our  two 
phases  as: 
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(Pi)°  _ 

100(p2)° 
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55 


(4.7) 


If  we  don't  want  to  make  the  assumption  that  the  solid  phase  is  always  in  equilibrium,  we  may  use 
relations  available  between  the  rate  of  change  of  concentration  of  the  solid  phase  and  the  relative 
humidity  of  the  gas  phase,  an  example  of  which  is  given  by  Norden  and  David9 . 


We  may  also  write  the  vapor  pressure-temperature  relation  for  the  vaporizing  (i  phase,  which 
Whitaker  gives  for  porous  media  as: 
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*  F- 

2<jSy  i 

,  ^K’ap 

i  n 

rpp/?i(r) 
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100  r.J 

(4.8) 


This  relation  gives  the  reduction  or  increase  in  vapor  pressure  from  a  curved  liquid  surface  resulting 
from  a  liquid  droplet  influenced  by  the  surface  interaction  between  the  solid  and  the  liquid,  usually 
in  a  very  small  capillary. 


In  many  cases,  the  Clausis-Clapeyron  equation  will  be  sufficiently  accurate  for  the  vaporizing 
species,  and  the  gas  phase  vapor  pressure  may  be  found  from: 


(P\Y  =  Pi  exp 
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&Kap 
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[  *1 

WO  rjjj 

(4.9) 


This  vapor  pressure-temperature  relation  is  only  good  if  we  have  the  lr  (uid  phase  present  in  the 
averaging  volume.  We  may  also  have  the  situation  where  we  only  have  the  solid  phase,  containing 
adsorbed  water,  and  the  gas  phase.  To  get  the  vapor  pressure  in  the  gas  phase  in  this  situation,  we 
will  use  the  sorption  isotherm,  and  assume  that  the  gas  phase  is  in  equilibrium  with  the  sorbed  water 
content  of  the  solid  phase. 
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We  can  use  any  isotherm  relation  where  we  have  the  solid  water  concentration  as  a  function  of 
relative  humidity.  The  equation  given  previously  is  one  example: 
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(4.10) 


From  a  curve-fit,  or  some  kind  of  equation  solver,  we  find  that  the  partial  vapor  pressure  can  be 
given  as  a  function  of  the  volume  fraction  of  water  in  the  solid  phase: 

(PiY  -  f{Ps’Pi’Ps’£ai)  a*  the  temperature  (T) ,  only  e0/  is  unknown.  (4.11) 
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5.  Mass  Transport  in  the  Gas  Phase 


The  volume  average  form  of  the  gas  phase  continuity  equation  was  found  to  be: 
|(eY(pT)1')+V-({pY)T(!>T))+^J/lifPy(!'y-iS>)-Hyl)d4+ij;V|PY(i>r-H>,)'rv<i4  =  0  (5.D 
and  the  species  continuity  equation  was  given  as: 

/  v 

^(eY<Pt)Y )+ V  *  ((Pi >T  )) + ^7  Pi  (v2  —  vi>)  •  =  v  •  ^  (py)  ©V  .  (5.2) 

where  the  dispersion  and  source  terms  were  dropped  from  the  equation. 

If  we  use  the  definition  of  the  mass  flux  from  one  phase  to  another  as: 

=  (5.3) 

or 

<m;v>  =  Py(pt  -  ») -V"  .  (5.4) 

with  the  same  form  for  the  mass  flux  from  the  solid  to  the  gas  phase,  the  gas  phase  continuity 
equation  may  be  rewritten  as: 


|-(Er{p,)1')+v{<pT)Y(vr))  =  (mk,)+{m„) 


For  the  two  species  (1— water,  and  2— air),  the  species  continuity  equations  are  written  (again 
dropping  the  source  and  dispersion  terms)  as: 


|(EY(pi)^)+V.((pi>^(vY))-K)-<<>-v--(p7)Y^ 


;(eY(P2>Y)+V-((p2>7(vY))  =  v',(Py)y2?v  ^ 


<p2>T 
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If  we  again  ignore  the  effects  of  the  dispersion  terms  in  the  diffusion  equations  derived  by 
Whitaker,  we  may  incorporate  an  effective  diffusivity  into  the  species  continuity  equations,  which 
are  now  given  as: 


|:Mpi)T)+vY(Pi)%)H'"k}“('"<") 


s.  \=v- 


'ieHL' 


|K(p2>T)+V.((p2^vr))=V.| 


(pr)T^v 


rM l 
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(5.8) 


(5.9) 


The  effective  diffusivity  will  be  some  kind  of  function  of  the  gas  phase  volume  e;  as  the  solid 
volume  and  the  liquid  volume  fractions  increase,  there  will  be  less  space  available  in  the  gas  phase 
for  the  diffusion  to  take  place.  We  might  try  to  define  the  effective  diffusivity  as: 


Deff  ~ 


A2^y  _ DaEy 


(5.10) 


X  X 

where  the  effective  diffusivity  Dtff  is  related  to  the  diffusion  coefficient  of  water  vapor  in  air  (Dn 
or  D  )  divided  by  the  effective  tortuosity  factor  x. 


An  example  of  a  good  relation  for  the  binary  diffusion  coefficient  of  water  vapor  in  air  is  given  by 
Stanish9  as: 


A2  - 


2.23 

f  T  V 
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U73.15J 

1.75 
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(5.11) 


To  simplify  matters,  one  could  assume  the  tortuosity  factor  to  be  constant,  and  let  the  variation  in 
the  gas  phase  volume  take  care  of  the  change  in  the  effective  diffusion  coefficient  as  the  volume 
available  for  gas  phase  diffusion  changes  with  solid  swelling  and/or  liquid  volume. 


Another  simplification  is  to  only  account  for  the  water  vapor  movement,  so  the  continuity  equation 
would  become: 
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6.  Gas  Phase  Convective  Transport 


It  is  important  to  include  forced  convection  through  porous  media  since  this  can  be  an  important 
part  of  the  transport  process  of  mass  and  energy  through  porous  materials  with  high  air 
permeability. 

It  is  not  necessary  to  modify  any  of  Whitaker's  derivations  for  the  gas  phase,  and  if  we  neglect 
gravity,  we  may  write  the  gas  phase  velocity  as: 

where  the  permeability  tensor  K, is  a  transport  coefficient. 

There  arc  other  methods  to  obtain  an  estimate  of  the  convective  velocity  of  a  gas  flow  through  a 
porous  material.  It  may  be  desirable  to  use  one  of  these  other  relations  to  obtain  the  volume 
average  form  of  the  gas  velocity. 


For  example,  we  could  start  directly  with  Darcy's  law: 
JLn 


VP  +  ~Vy  =0 
K  y 


(6.2) 


and  assume  that  for  the  dry  porous  material  we  have  available  the  experimental  measurement  of  the 
specific  permeability  coefficient  K,  and  then  modify  it  to  account  for  the  decrease  in  gas  phase 
volume  as  the  solid  phase  swells  and/or  the  liquid  phase  accumulates.  We  could  make  the  variation 
in  K  a  linear  function  of  the  gas  phase  volume,  which  has  been  an  approach  used  by  Stanish  et.  al.u 

1  (6.3) 


£fdry 


This  is  a  very  simple  model,  and  may  be  improved  upon.  In  the  book  by  Dullien12,  there  are  a 
variety  of  relationships  for  how  K  varies  with  porosity;  some  of  those  relations  may  be  more 
realistic  for  our  purposes.  We  could  also  relate  the  change  in  the  permeability  to  the  effective 
tortuosity  function  t,  which  also  has  the  same  factors  related  to  the  decrease  in  gas  phase  volume, 
and  change  in  geometry,  that  we  need  to  account  for  the  Darcy's  law  relation  for  convective  gas 
flow. 
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7.  Liquid  Phase  Convective  Transport 


Whitaker's  derivation  for  the  convection  transport  of  the  liquid  phase  is  the  one  of  the  most 
complicated  parts  of  his  general  theory.  He  accounts  for  the  capillary  liquid  transport,  which  is 
greatly  influenced  by  the  geometry  of  the  solid  phase,  and  the  changeover  from  a  continuous  to  a 
discontinuous  liquid  phase.  His  eventual  transport  equation,  which  gives  an  expression  for  the 
liquid  phase  average  velocity  is  quite  complicated,  and  depends  on  several  hard-to-obtain  transport 
coeffients.  The  final  equation  is  given  as: 


•  [^e^Ep  +  k(T)V(T)  -  (pp  -  py )] 


(symbol  definitions  given  in  nomenclature  table) 


(7.1) 


One  advantage  of  Whitaker's  derivation  is  that  it  is  almost  completely  independent  of  the  other 
transport  equation  derivations.  This  should  mean  that  we  may  use  another  expression  for  the  liquid 
phase  velocity  if  we  find  one  that  is  more  amenable  to  experimental  measurement  and  verification. 


An  example  of  an  equation  which  is  more  empirical  is  again  given  by  Stanish11.  The  velocity  is 
assumed  proportional  to  the  gradient  in  pressure  within  the  liquid.  The  pressure  in  the  liquid  phase 
is  assumed  to  be  the  sum  of  the  gas  pressure  within  the  averaging  volume  minus  the  capillary 
pressure  (P): 


If  we  use  a  relation  of  this  kind  it  is  necessary  to  obtain  an  equation  for  the  capillary  pressure  as  a 
function  of  the  fraction  of  non-solid  volume  occupied  by  the  liquid  phase,  as  well  as  a  relation  for 
the  variation  in  the  permeability  coefficient  as  a  function  of  liquid  phase  volume  fraction.  It  is  also 
necessary  to  to  determine  when  the  liquid  phase  becomes  discontinuous  so  that  liquid  flow  ceases  at 
that  point.  These  types  of  relations  can  be  determined  experimentally  for  materials  of  interest,  or 
they  may  be  found  in  the  literature  for  quite  a  wide  variety  of  materials. 
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8.  Summary  of  Modified  Transport  Equations 

The  set  of  modified  equations  which  describe  the  coupled  transfer  of  heat  and  mass  through 
hygroscopic  porous  materials  are  summarized  below. 


Total  thermal  energy  equation  (8.1) 

M 

<p)C,$22+  +Pp(c,)p(%)  •V(r>+A*vap(mfc)+0(mt,)+(a  +  AA,»PK'"*v>  =  vW-v<r>) 

.  1=1 


Liquid  phase  equation  of  motion 
(^p)  =  -  ^  V(<Pi)Y+(P2)Y-^) 


Liquid  phase  continuity  equation 

if  +  V '(%)  +  ^  (%'*)' + ^  (% 

which  may  be  rewritten  as: 


-w2)-iipotM  =  0 


«m/v  )-(<)) 

Pp 


=  0 


Gas  phase  equation  of  motion 

(^)  =  -f^LV(<A)Y+{P2>Y) 


Gas  phase  continuity  equation 

|(er(pr)T)^{{p *)*(«,))-<*>+<*.) 


(8.2) 

(8.3) 

(8.40 

(8.5) 

(8.6) 


39 


Gas  phase  diffusion  equations 


Solid  phase  density  relations 

w-w+w 

Pi  =  ec lPl 
p2  =  (1-E0l)P5 
eaS  +  £aL  “1 

Solid  phase  continuity  equation 

(£e  (Per )° ) + ^  •  ((p0  )°  ( vff  >) + (msl ) + (m^ )  =  0 

Solid  phase  equation  of  motion  (for  one  dimensional  geometry) 

Solid  phase  diffusion  equation  (for  vaporizing  component) 

{  y 

|:(£o(pi)0)+V-((p1}0(v1))+(mJ,)  +  {<)  =  V-'(p0)a©0V  ^5-  * 
Volume  constraint 

eo(0+ep(f)+eY(0-l 

Thermodynamic  relations 

<Pi>Y  =(Pi)Y/?i<r> 

<p2)Y=(p2)Y/?2<r) 

(Pr)Y=<Pi>Y+<P2>Y 

(Pyf  a{PlY +Wit 


(8-9) 

(8.10) 

(8.11) 

(8-12) 


(8.13) 


(8.14) 


(8.15) 


(8.16) 


(8.17) 

(8.18) 

(8.19) 

(8.20) 


40 


If  any  liquid  phase  is  present,  vapor  pressure  is  given  by: 
(PiY  =  Pi  exp 


2^  > 

,  ^Kap 

i  n 

/P^iiT); 

[(r)'j;J 

or 


{PiY  =  Pi  expj- 


A 


>ap 


_1 _ 1_ 

(T)  Z 


(8.21) 


If  the  liquid  phase  is  not  present,  and  the  liquid  component  is  desorbing  from  the  solid,  the  reduced 
vapor  pressure  in  equilibrium  with  the  solid  must  be  used.  This  relation  may  be  determined  directly 
from  the  sorption  isotherm  for  the  solid: 

{PiY  =  f(PsiPi’Ps’EoL)  at  temperature  { T )  ,  only  zaL  is  unknown  .  (8.22) 


Sorption  relations  (volume  average  solid  equilibrium) 

(  ^ 


Qsi  (J/ kg)  =  0.195 


1- 


M 


0.2  + 


iEl 


1.05 


Ps 


(8.23) 


(P2)°  (^”eo/)Pi 

0.55^- 

l  J 

fo.25+^] 

•+■ 

.1  A  J 

l  A  J 

(8.24) 


This  is  a  total  of  20  main  equations  and  20  unknown  variables,  which  should  allow  for  the  solution 
of  the  set  of  equations  using  numerical  methods.  The  20  unknown  variables  are: 


si )» (^jv  )>  (™iv  )*  Qsl 


(p,)T.{ft)T.{p2)7.(pr>Y.<p1)T.<p2)r 

(py)°.(p,)”.<p2><’ 
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9.  Simplification  of  Total  Equation  Set 

The  total  equation  set  is  quite  complicated.  There  are  several  ways  to  simplify  the  set  of  equations, 
yet  still  include  the  modes  of  energy  and  mass  transport  which  are  important  in  various  situations. 
First,  we  will  include  the  the  gas  phase  convection  and  the  liquid  phase  convection  terms,  but  make 
several  assumptions  about  the  diffusion  processes  that  will  simplify  the  set  of  equations 
significantly.  Then  we  will  further  simplify  the  set  of  equations  to  examine  the  case  where  gas 
phase  convection  is  not  present,  but  liquid  capillary  transport  is  important.  Finally,  we  will  simplify 
the  set  of  equations  to  the  situation  where  mass  is  only  tranported  through  diffusion  in  the  gas 
phase. 

Simplified  equations  which  include  liquid  and  gas  phase  convection 
The  total  drying  equation  set  will  be  simplified  by  making  several  assumptions: 

1)  We  will  only  use  the  continuity  equations  to  account  for  the  mass  transfer  of  the  component 
which  crosses  phase  boundaries  (water).  This  means  we  will  only  use  the  species  continuity  and 
diffusion  equations  for  component  1. 

2)  We  will  use  the  diffusion  coefficients  for  water  in  the  solid  and  water  in  air,  and  ignore  the 
counterdiffusion  of  air  through  the  water  vapor,  etc.  We  will  assume  these  diffusion  coefficients 
are  constant. 

3)  For  the  diffusion  equations,  we  will  assume  that  the  density  of  the  inert  phase  (solid  or  air)  is 
constant  during  the  diffusion  process. 

4)  We  will  ignore  the  swelling  velocity  of  the  solid  phase,  but  will  the  include  the  change  in 
volume  of  the  solid  as  water  is  absorbed, 

5)  We  will  assume  transport  only  in  one  dimension  (x). 


Total  thermal  energy  equation 


PeMpW 

W.+ 

+Qt(”lsi) 

__3_ 

r,  wi 

+(<>), (Pi»l) 

3x 

+(Qf  +  AAVfl/,){mJV)_ 

9x 

ly&j 

(9.1) 


Liquid,  phase  equation  of  motion 


42 


Liquid  phase  continuity  equation 

l»(atfehM=0 

dt  3jc'  p'  Pp 

Gas  phase  equation  of  motion 

Gas  phase  continuity  equation 

|(eT(pY)Y)+£(<Pr)1,{v1»))=K>+{/(i„} 


Gas  phase  diffusion  equation 


|:(ey<Pi)Y)+|;((pi  }Y{v)H<M^)  =  a; 


*  ((pl>Y) 


(9.3) 


(9.4) 


(9.5) 


(9.6) 


Solid  phase  density  relations 
(Po)<T={Pi)<,+(P2)° 

Pi  =  EoZ.Pl 
p2  =  (1  -e0l)Ps 

EoS  +  eoL  = 1 

Solid  phase  continuity  equation 
Jt(eo(PoT)+(^si)+{<)  =  0 


Solid  phase  diffusion  equation 


Volume  constraint 

eo(/)  +  ep(r)+eY(0  =  l 


(9.7) 

(9.8) 

(9.9) 

(9.10) 

(9.11) 


(9.12) 

(9.13) 
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Thermodynamic  relations 

(pi)r=<Pi>r  *i(r) 

(P2>r=(P2>’,«2(r> 

faY-W+W 

(PY)r=<Pl)T+<P2>T 


If  any  liquid  phase  is  present,  vapor  pressure  is  given  by: 


(9.14) 

(9.15) 

(9.16) 

(9.17) 


(9.18) 

(9.19) 


If  the  liquid  phase  is  not  present,  and  the  liquid  component  is  desorbing  from  the  solid,  the  reduced 
vapor  pressure  in  equilibrium  with  the  solid  must  be  used.  This  relation  may  be  determined  directly 
from  the  soiption  isotherm  for  the  solid: 

( PiY  =  /(Pj-P/>Pj>eo/)  at  the  temperature  (T) ,  only  £o/  is  unknown  (9.20) 


Sorption  relations  (solid  equilibrium) 

(  \ 


Qi  (J/ kg)  =  0.195 


1- 


<p.n 

1 

1 

1 

p‘  J 

[o.2+M] 

"T" 

A  Ps  > 

l  Ps  J 

(9.21) 


(Pl)°  =  £q/P/ 
(P2)0  (^-e0/)Pj 


=  R< 


0.55 


Ps 


y\ 


(o.25+M1 

f,. 25-^1 

l  Ps  J 

l  Ps  J 

(9.22) 
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Transport  Coefficients  and  Mixture  Properties 


Kff  -  eo 


in  '  "in  +ct +e  an.:  rv: 
P1+P2  )  ^  y{  P1+P2 


(9.23) 


%=■ 


<P>  =  e0(<Pi)°  • +{p2)c)+e0(pp)|i + eY({pi)Y +(p2)Y) 

„  «a[<Pi>a  (cpX  +<P2>°  (^>2] +epPp(«pXi + «r[<ft>r  (">X  +<P*>Y  (cj-)2] 

<P> 


(9.24) 


(9.25) 


(9.26) 


This  simplified  equation  set  is  now  19  equations  for  19  variables,  which  is  enough  to  solve  the  set 
of  equations  for  the  following  unknown  variables: 


eo.£p.er(?p).(?y),CO. 

{  ^sl )  .("*«  )>  (f»lv  )  J  Qst 

(Py)\{PiY '(PlY  >(Py)\(pl)y  >(PlY , 

{PyY  >(Pl)y  ’(PlY 
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Simplified  equations  which  neglect  gas  phase  convection 


The  gas  phase  convection  is  now  neglected,  so  that  the  air  within  the  pore  spaces  is  assumed  to  be 
stagnant,  and  mass  and  energy  are  transported  within  the  gas  phase  only  by  diffusion. 


Total  thermal  energy  equation 


(P  >C, 


m 

d, 


+ 


[pf(c/')s(v 


)]f+l 


^ap(K) 

+Qi{/”si) 

+fe  +  Artvap){m,rv) 


d_ 

dx 


<eff 


dx 


(9.27) 


Liquid  phase  equation  of  motion 


(9.28) 


Liquid  phase  continuity  equation 

3,  '  (K)-K»- 

dt  atcvP/  Pp 

Gas  phase  continuity  equation 
|(8y(Py)Y)  =  {^v)+R) 

Gas  phase  diffusion  equation 

|(ey(p,y)-<»„)-K)=|[%|({Ply)l 


(9.29) 


(9.30) 


(9.31) 


Solid  phase  density  relations 

(po)°  ={Pi)0+{P2)C 

Pi  = 6olPl 

p2  =(l-eo£,)ps 

eoS+eaL  =1 


Solid  phase  continuity  equation 
|:(£a(po)CI)-t-(^)+{^)  =  0 


Solid  phase  diffusion  equation 

|-(£cT{pi)a)+{^>+<^)  =  -^ 


©a 


(9.32) 

(9.33) 

(9.34) 

(9.35) 

(9.36) 


(9.37) 
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Volume  constraint 

ea(0+ep(0+eY(f)  =  l 
Thermodynamic  relations 

Wr=<p.>T*i<n 
(pj>,=(p2>T^(n 
(p,)1'=(p.>,+(p2>T 
(Pt)Y  =  (Pi>Y+<P2>1' 


If  any  liquid  phase  is  present,  vapor  pressure  is  given  by: 


{PiY  =  Pi  exP 


y  0  _  \ 

2cPy 

l/Pp^iOO, 

AAy, 


'2L 


JL  JL 

(T)~Z 


or  {PiY  =  Pi°  exp 


’A^vap 

r  i  nil 

L  *1 

loo  rJJ 

(9.38) 

(9.39) 

(9.40) 

(9.41) 

(9.42) 


(9.43) 

(9.44) 


If  the  liquid  phase  is  not  present,  and  the  liquid  component  is  desorbing  from  the  solid,  the  reduced 
vapor  pressure  in  equilibrium  with  the  solid  must  be  used.  This  relation  may  be  determined  directly 
from  the  sorption  isotherm  for  the  solid: 

{PiY  =  fiPs’PbPs >8o/)  at  temperature  (r)  ,  only  e0/  is  unknown  (9.45) 


Sorption  relations  (solid  equilibrium) 

r 


Qi  (J/kg)  =0.195 


Ps 


).2+^ 

Pi 


r  w 


1.05- 


V 


Ps 


(9.46) 


<Pi)q  =  eq<P/ 
(p2)°  (1  —  Ea/)Pi 


=  R 


'S 


0.55 


WY) 

1 

4_ 

1 

A  J 

fo.25+M] 

i 

[L25»T] 

.1  Ps  J 

l  Ps  J 

(9.47) 
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Transport  Coefficients  and  Mixture  Properties 


Kff  —  £<j 


P1  +  P2 


‘£p^p 


+  e. 


^lPl  +  ^2p2 

Pl+P2  . 


^  ®«eY 

**—?■ 


(p)  ■  eo«Pi)° + (p2>°)+e3(pp)p + ey«p1>Y+<P2>T) 

gCT[<Pi)°(gJ.)1+<P2>°(cP)2]+gpPp(cp)B-»-eT[<p1)T(cp)l+(p2)T(cf,)2] 


c„  = 


<P> 


(9.48) 

(9.49) 

(9.50) 

(9.51) 


This  simplified  equation  set  is  now  18  equations  for  18  variables,  which  is  enough  to  solve  the  set 
of  equations  for  the  following  unknown  variables: 


So»Ep»£y»(^p)»(7,)» 
{fist ),  (jflyp),  (miv),  Qyl 


{pyf  >{PlY>{P2y>{  Py)\{Pl)\{P2)\ 

{Pr)Y.(Pi)Y.(P2>Y 
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Simplified  equations  which  neglect  both  liquid  and  gas  phase  convection 


Both  liquid  phase  capillary  transport,  and  gas  phase  convection  are  neglected,  so  that  mass  and 
energy  are  transported  within  the  gas  phase  only  by  diffusion,  and  the  liquid  which  condenses  at  a 
particular  point  within  the  material  does  not  wick  or  flow  away  from  that  point,  but  accumulates 
over  time.  The  liquid  may  be  absorbed  into  the  solid  phase,  or  it  may  evaporate  at  a  later  time,  but 
it  may  only  be  transported  through  the  structure  by  going  into  the  vapor  phase. 


Total  thermal  energy  equation 


<p>c^+i 


^Kapi^tv) 

+Qi{^sl) 

+  +  A/ \0p  ){ ^ 


Liquid  phase  continuity  equation 

3eP  ,  «*%)-(*»)) .  Q 

7)t  pp 

Gas  phase  continuity  equation 
— (ay(pY)  j  =  {»J/v)  +  (wJV) 

Gas  phase  diffusion  equation 


J~(ey(Pl)Y)  "{'"/v)  -{'O  = 


(9.52) 


(9.53) 


(9.54) 


(9.55) 


Solid  phase  density  relations 

W-W+fey 

Pi  =  eaZ,PZ. 

P2  =0“eaZ.)Ps 

EaS' +  eaL  =  1 


Solid  phase  continuity  equation 
^(Eo{Po>a)+<^/)+<'«,v>=0 

Solid  phase  diffusion  equation 


o(pi  )0)+(w,/) +(/»*,) 


(9.56) 

(9.57) 

(9.58) 

(9.59) 

(9.60) 


(9.61) 
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Volume  constraint 

eo(0+£p(r)+£Y(0=l 

(9.62) 

Thermodynamic  relations 

<p.)H  Pi)T«i(r) 

(9.63) 

(9.64) 

<Pt)T  =<Pi>T+(Pj)T 

(9.65) 

(pt)T={a),+(P2>1' 

(9.66) 

If  any  liquid  phase  is  present,  vapor  pressure  is  given  by: 


(PiY  =  Pi  exp 
or  (pi)y  =  p\  exp 


7  "\ 

2oW 

A h^(  1  lx 


l  (T)  Te 


1  nil 

. 

l(r>  r.JJj 

(9.67) 

(9.68) 


If  the  liquid  phase  is  not  present,  and  the  liquid  component  is  desorbing  from  the  solid,  the  reduced 
vapor  pressure  in  equilibrium  with  the  solid  must  be  used.  This  relation  may  be  determined  directly 
from  the  sorption  isotherm  for  the  solid: 

(PiY  =  f{Ps>PhPs'eat)  at  the  temperature  (T) ,  only  e0/  is  unknown  (9.69) 


Qi  (J/kg)  =  0.195 


't.M" 

Ps 


1 


0.2  + 


(jI 


1.05 


_i£ll 

Ps 


(9.70) 


(Pi)°  _  ?o/P/  -R 

f0  55W) 

1  1 

X 

(p2)°  f 

l  A  J 

fo.25+^] 

A  a  j 

T 

f,25-^l 

l  ^  , 

(9.71) 
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Transport  Coefficients  and  Mixture  Properties 


Keff 


=  e  (!^2i±h£i 


%=■ 


V  Pl+P2 


£p/^  +  £y 


&1P1  +  kqffi 

v.  Pl+P2  , 


(p)  =  ea(<Pi)0+(p2)°)+ea(pp)P  +  £Y(<Pi>Y+(p2>Y) 


c„  = 


£q[<Pl)q  {CP  )i  +  <p2>°  {cp  )2]  +  £ppp  (cp  )g  +  ey  [<Pl)Y(cp  )t  +  <p2  )Y  (<7>  )2 


<P) 


(9.72) 

(9.73) 

(9.74) 

(9.75) 


This  simplified  equation  set  is  now  17  equations  for  17  variables  (the  sorption  equation  and  the 
solid-vapor  pressure  equation  are  a  duplicate),  which  are  enough  to  solve  the  set  of  equations  for 
the  following  unknown  variables: 


£or.£p»£y»\r /, 

( ^hl )» ( )  >  { Mfa  )  ’  Qsl 


(p7)r.{A>T.<ft)T.(pT>Y.(p,>T.(p2)r, 

(Py)T.(Pi}1'.{P2>T 
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10.  Comparison  with  Previously-Derived  Equations 

The  simplified  sytems  of  partial  differential  equations  given  in  the  previous  chapter  still  contain 
many  equadons  with  a  large  number  of  unknown  variables.  Even  for  the  simplified  case  of  vapor 
diffusion,  the  system  of  equations  is  quite  confusing,  and  it  is  difficult  to  verify  their  accuracy, 
other  than  by  checking  for  dimensional  consistency.  One  way  of  checking  their  validity  is  to  see  if 
they  simplify  down  to  more  well-known  diffusion  equadons  for  the  transport  of  water  vapor  in  air 
through  a  porous  hygroscopic  solid.  Such  a  system  of  equadons  has  been  well  documented  by 
Henry12,  Norden  and  David9 ,  and  Li  and  Holcombe13 ,  who  have  used  them  to  describe  the 
diffusion  of  water  vapor  through  a  hygroscopic  porous  material. 

We  will  make  the  same  assumptions  used  by  these  previous  workers,  and  attempt  to  transform  the 
system  of  equations  for  the  case  of  vapor  diffusion  (no  liquid  or  gas  phase  convection)  to  their 
system  of  equations.  For  completeness,  we  will  also  need  to  write  the  various  equations  in  terms  of 
the  variables  and  units  used  their  work. 


The  major  simplifying  assumptions  are:  1)  there  is  no  liquid  or  gas  phase  convection,  2)  there  is  no 
liquid  phase  present,  3)  the  heat  capacity  of  the  gas  phase  can  be  neglected,  4)  the  volume  of  the 
solid  remains  constant  and  does  not  swell,  5)  the  solid  and  gas  phase  volume  fractions  are  both 
constant,  6)  the  thermal  conductivity  tensor  may  be  expressed  as  a  constant  scalar  thermal 
conductivity  coefficient,  7)  the  gas  phase  diffusion  coefficient  is  constant,  7)  the  transport  is  one¬ 
dimensional  (x-direction). 

The  total  thermal  energy  equation  becomes: 

(pK  ^-+( a + a vX“»> = *-K  ■' m)  ao.D 

or 

(p)cp  ^+(C/+ A Kap  )(<-  >  =  Kff  pp 

The  gas  phase  continuity  equation  becomes: 

‘rf«Pr»K> 


(10.2) 


(10.3) 


The  gas  phase  diffusion  equation  (component  1-water  vapor): 


eyf((pl>Y)"(^) 


32(Pi)Y 

Ox2 


(10.4) 


(10.5) 
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The  solid  phase  continuity  equation  (component  1 -water): 

e»^((pi)<’)+{m„)  =  °  (10.6) 


For  the  solid  phase  diffusion  equation  (component  1 -water)  we  assume  that  the  diffusional 
transport  through  the  solid  phase  is  insignificant  compared  to  the  diffusion  through  the  gas  phase,  so 
the  diffusion  equation  reduces  to  the  continuity  equation: 


e0-|((pl>°)+(^>  =  V' 


<Pc>£ 


Pi 


,<Po>° 


=  0 


(10.7) 


Volume  fraction  constraint 

£y  +  £cf  =  1  £(J  =  1  —  £y 

Thermodynamic  relations 

<Pi)Mp.)Y  *.<!•> 

{piY  Ri(T) 
(Py)T=<Pi)Y  +  <P2>Y 
(pyf  =(PlY +(P2)y 


(10.8) 

(10.9) 

(10.10) 
(10.11) 
(10.12) 


We  made  the  assumption  that  the  mass  transport  through  the  solid  phase  is  negligible  compared  to 
mass  transport  through  the  gas  phase.  This  is  reasonable  since  the  diffusion  coefficient  for  water  in 
a  solid  is  always  much  less  than  the  diffusion  coefficient  of  water  vapor  through  air.  We  thus  only 
have  accumulation  of  water  in  the  solid,  and  the  solid  acts  as  a  source  or  sink  for  water  vapor. 


We  can  combine  the  continuity  equations  for  water  (component  1)  for  both  phases  by  connecting 
the  phase  equations  through  the  mass  flux  from  the  solid  to  the  gas  phase: 


ea|((Pi>°)  +  K)  =  0 
er|((pi>Y)-«>  =  %^r- 


(10.13) 

(10.14) 


Bt 


(<Pi  >Y)-(^> 


which  we  may  rewrite  in  terms  of  the  gas  phase  volume  fraction  as: 


('-ev)|((P>)°)  +  eY|((P>)Y)=®^ 


i2M 

3x2 


(10.15) 


(10.16) 
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Through  these  various  assumptions,  we  have  reduced  our  large  equation  set  down  to  two  main 
equations  for  the  energy  balance  and  the  mass  balance: 


(p)cP  +  (fir +  Mvap  )(4v)  =  Kff 

(1-Er)|((Pl)°)+e4«pi>Y)=%¥ 


(10.17) 

(10.18) 


To  make  the  comparison  with  the  existing  equations  of  Henry12,  Norden  and  David9 ,  and  Li  and 
Holcombe13, easier,  we  can  rewrite  the  intrinsic  phase  averages  in  terms  of  the  concentration  of 
water  in  the  solid  ( CF )  and  concentration  of  water  in  the  gas  phase  (C): 


_  mass  of  water  in  solid  phase  _  /nlp  _ 

F  solid  phase  volume  Va  ° 


(10.19) 


C  mass  of  water  in  gas  phase  miy 


gas  phase  volume 


(10.20) 


Since  the  definition  of  intrinsic  phase  average  gives  the  same  quantity  as  the  true  point  value,  we 
may  use  the  fact  that 

W-toF-cv  (1°'21> 

(Pl>M  cf=c  (10-22> 

to  rewrite  the  mass  balance  equation  as: 

ar-^-S— <0  (,0-23) 


(l-e7)^+ev~  =  2>. 


We  can  rewrite  the  effective  diffusion  coefficient  by  using  the  diffusion  coefficient  for  water  vapor 
in  air  modified  by  the  gas  volume  fraction  and  the  tortuosity  of  the  gas  volume  fraction: 


(i 


\  dCF  dc  Daey  a2c 

Ey  I  ^  “H  Gy  ~  “ 

v  dt  7  8/ 


(10.24) 


The  thermal  energy  equation 

<P  )CP  ^ + (Qt +  )R )  =  Kff 

may  also  be  modified  by  recognizing  that  the  mass  flux  term  is  contained  in  the  solid  phase 
continuity  equation: 

e°  Jf  ((pl>° ) + <^w> =  °  ^ 

so  that  the  thermal  energy  equation  may  be  rewritten  as: 

8CV  ,  a2(r) 


(10.25) 


(10.26) 


(p)cp  ^p-(a = v 


(10.27) 
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If  we  go  back  to  our  definitions  for  the  mass  fraction  weighted  average  heat  capacity, 

c  -  _ i! _ 

(P> 

and  spadal  average  density, 

<p)=e„f{p,}XX<p.)T 

;=1  M 

the  thermal  energy  equation  may  be  rewritten  as: 

e,[<Pi)'’(c,)1+(P2}0(cp)2]+eY[(p1)1,(cp)1+(p2>T(cp)2]  ^ 

-(a+A 


(10.28) 


(10.29) 


(10.30) 


(10.31) 


If  we  make  the  assumption  that  the  heat  capacity  of  the  gas  phase  is  negligible,  then  the  thermal 
energy  equation  becomes: 


ea[<Pi)°(cp)1+<P2)0(^)2]  +  =  * 

or  dividing  through  by  the  solid  volume  fraction: 

[(Pi)°(^)1+<P2><’(c,)2]^— = 

For  consistent  nomenclature  with  Li  and  Holcombe14  we  will  write  the  effective  thermal 
conductivity  keff  as  K. 

We  can  also  define  a  volumetric  heat  capacity  Cv  as: 


(10.31) 


(10.32) 


UnitSfor(py)°(Cp). 


(10.33) 


The  final  thermal  energy  equation  reduces  to: 

c.fi 


(10.34) 
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(10.35) 


The  two  simplified  equations  for  the  mass  and  energy  balance  are  thus: 

i  scr  . .  sc  djl y  a2c 

-e^-3r+e’¥=_ra?- 

cv^i-(a+'Uw)^t=« 


(10.36) 


These  two  simplified  equations  are  very  encouraging,  since  they  are  exactly  the  same  as  previous 
equations  derived  by  Henry12,  Norden  and  David9 ,  and  Li  and  Holcombe13  for  describing  the 
diffusion  of  water  vapor  through  a  hygroscopic  porous  material.  In  their  equations  they  define  the 
heat  of  sorption  from  the  vapor  phase  into  the  solid  (which  is  the  opposite  of  the  heat  of  desorption 
which  we  used)  as: 


lk  =  [Q,  +  Ahyap ) 

so  that  their  equations  are: 


y  32C 
X  dx2 


(10.37) 

(10.38) 


(10.39) 
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11.  Conclusions 


Whitaker's  theory  of  coupled  heat  and  mass  transfer  through  porous  media  was  modified  to  include 
hygroscopic  porous  materials  which  can  absorb  liquid  into  the  solid  matrix.  The  system  of 
equations  described  in  this  report  should  make  it  possible  to  evaluate  the  time-dependent  transport 
properties  of  hygroscopic  and  non-hygroscopic  clothing  materials  by  including  many  important 
factors  which  are  usually  ignored  in  the  analysis  of  heat  and  mass  transfer  through  textile  materials. 
The  set  of  equations  allows  for  the  unsteady  capillary  wicking  of  sweat  through  fabric  structure, 
condensation  and  evaporation  of  sweat  within  various  layers  of  the  clothing  system,  forced  gas 
phase  convection  through  the  porous  structure  of  a  textile  layer,  and  the  swelling  and  shrinkage  of 
fibers  and  yams  as  they  absorb/desorb  liquid  water  and  water  vapor. 

The  simplified  set  of  equations  for  heat  and  mass  transport,  where  mass  transport  occurs  due  to 
diffusion  within  the  air  spaces  of  the  porous  solid,  was  shown  to  reduce  to  the  well-known  coupled 
heat  and  mass  transfer  models  for  hygroscopic  fabrics,  as  exemplified  by  the  work  of  Li  and 
Holcombe14. 

Work  is  underway  to  develop  a  numerical  code  to  solve  the  various  sets  of  equations  by  standard 
numerical  methods  based  on  Patankar's  control  volume  approach14.  The  numerical  method  is 
similar  to  those  developed  by  Whitaker  and  Vafai15,  and  Tao,  Besant,  and  Rezkallah16’17,18,  which 
have  had  success  in  modeling  the  unsteady  coupled  heat  and  mass  transfer  process  in  fibrous 
insulation  materials. 
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